What does it mean to be correlated?
I created a simulation study to observe the co-evolutionary outcome of the newt-snake interaction with different genetic architectures (GAs) in a spatial setting. I hypothesized that we would see an interaction (co-evolutionary arms race) between newt and snake phenotype under some GA combinations when newts and snakes were evolving over geographical space. Each GA is paired with another GA creating 16 combinations.
Experiment values:
many of the experiments did not run to completion (i will run this part 3 times)
## All cor, lit, and grid files exist!
## This program will now end!
## All cor, lit, and grid files exist!
## This program will now end!
## All cor, lit, and grid files exist!
## This program will now end!
## All cor, lit, and grid files exist!
## This program will now end!
## Group.1 x
## 1 1e-08_0.005_1e-08_0.005 0.124241445
## 2 1e-08_0.005_1e-09_0.0158 -0.073054703
## 3 1e-08_0.005_1e-10_0.05 -0.028004664
## 4 1e-08_0.005_1e-11_0.158 -0.156932495
## 5 1e-09_0.0158_1e-08_0.005 -0.333876391
## 6 1e-09_0.0158_1e-09_0.0158 0.169018200
## 7 1e-09_0.0158_1e-10_0.05 0.043688199
## 8 1e-09_0.0158_1e-11_0.158 -0.205526939
## 9 1e-10_0.05_1e-08_0.005 -0.034107364
## 10 1e-10_0.05_1e-09_0.0158 -0.166204923
## 11 1e-10_0.05_1e-10_0.05 -0.003520529
## 12 1e-10_0.05_1e-11_0.158 -0.243898057
## 13 1e-11_0.158_1e-08_0.005 -0.067010232
## 14 1e-11_0.158_1e-09_0.0158 -0.255090902
## 15 1e-11_0.158_1e-10_0.05 0.108156482
## 16 1e-11_0.158_1e-11_0.158 -0.034887799
## Group.1 x
## 1 1e-08_0.0158_1e-08_0.0158 -0.0009979573
## 2 1e-08_0.0158_1e-09_0.05 -0.3831566904
## 3 1e-08_0.0158_1e-10_0.158 0.1964659133
## 4 1e-08_0.0158_1e-11_0.5 -0.5891484283
## 5 1e-09_0.05_1e-08_0.0158 -0.9085769020
## 6 1e-09_0.05_1e-09_0.05 0.0100866952
## 7 1e-09_0.05_1e-10_0.158 -0.0468457273
## 8 1e-09_0.05_1e-11_0.5 -0.2917227943
## 9 1e-10_0.158_1e-08_0.0158 -0.6769856050
## 10 1e-10_0.158_1e-09_0.05 -1.1172741054
## 11 1e-10_0.158_1e-10_0.158 0.1573566841
## 12 1e-10_0.158_1e-11_0.5 -0.2728443554
## 13 1e-11_0.5_1e-08_0.0158 -0.3926328023
## 14 1e-11_0.5_1e-09_0.05 -1.4163172392
## 15 1e-11_0.5_1e-10_0.158 -0.8864611875
## 16 1e-11_0.5_1e-11_0.5 -0.1341703587
## Group.1 x
## 1 1e-08_0.05_1e-08_0.05 -0.4040913
## 2 1e-08_0.05_1e-09_0.158 -0.5262725
## 3 1e-08_0.05_1e-10_0.5 -0.5856175
## 4 1e-08_0.05_1e-11_1.58 -0.6708966
## 5 1e-09_0.158_1e-08_0.05 -0.4935271
## 6 1e-09_0.158_1e-09_0.158 0.3489808
## 7 1e-09_0.158_1e-10_0.5 -0.5490028
## 8 1e-09_0.158_1e-11_1.58 -0.5515514
## 9 1e-10_0.5_1e-08_0.05 -0.3452146
## 10 1e-10_0.5_1e-09_0.158 -0.5391440
## 11 1e-10_0.5_1e-10_0.5 -0.5478391
## 12 1e-10_0.5_1e-11_1.58 -0.4274895
## 13 1e-11_1.58_1e-08_0.05 -0.2236880
## 14 1e-11_1.58_1e-09_0.158 -0.2568742
## 15 1e-11_1.58_1e-10_0.5 -0.3409935
## 16 1e-11_1.58_1e-11_1.58 -0.3736108
## Group.1 x
## 1 1e-08_0.158_1e-08_0.158 -0.63330968
## 2 1e-08_0.158_1e-09_0.5 -0.50253898
## 3 1e-08_0.158_1e-10_1.58 -0.85901791
## 4 1e-08_0.158_1e-11_5 -0.09325217
## 5 1e-09_0.5_1e-08_0.158 -0.05162931
## 6 1e-09_0.5_1e-09_0.5 -0.21660078
## 7 1e-09_0.5_1e-10_1.58 -0.21723654
## 8 1e-09_0.5_1e-11_5 -0.87978708
## 9 1e-10_1.58_1e-08_0.158 -0.67831601
## 10 1e-10_1.58_1e-09_0.5 -0.10270499
## 11 1e-10_1.58_1e-10_1.58 -0.04610612
## 12 1e-10_1.58_1e-11_5 -0.64959493
## 13 1e-11_5_1e-08_0.158 -0.14630367
## 14 1e-11_5_1e-09_0.5 0.52834103
## 15 1e-11_5_1e-10_1.58 0.11150282
## 16 1e-11_5_1e-11_5 -0.58431132
## [1] "pattern 1e-08_0.005_1e-11_0.158_1"
## [1] "Cor between average snake pheno and local cor 0.853760572778879"
## [1] "Cor between average newt pheno and local cor 0.362671403041944"
## [1] "Cor between average dif pheno and local cor 0.879134338226174"
## [1] "Cor between newt pheno and snake 0.625012008488808"
## [1] "pattern 1e-08_0.005_1e-10_0.05_2"
## [1] "Cor between average snake pheno and local cor 0.77674759413511"
## [1] "Cor between average newt pheno and local cor 0.814466434187287"
## [1] "Cor between average dif pheno and local cor 0.235960198145116"
## [1] "Cor between newt pheno and snake 0.919844692645017"
## [1] "pattern 1e-11_0.158_1e-08_0.005_2"
## [1] "Cor between average snake pheno and local cor 0.376721378075875"
## [1] "Cor between average newt pheno and local cor 0.204801959367532"
## [1] "Cor between average dif pheno and local cor -0.0723515773520942"
## [1] "Cor between newt pheno and snake 0.902384185258089"
## [1] "pattern 1e-08_0.0158_1e-08_0.0158_2"
## [1] "Cor between average snake pheno and local cor -0.168012931704619"
## [1] "Cor between average newt pheno and local cor -0.227831237202084"
## [1] "Cor between average dif pheno and local cor 0.321188254748969"
## [1] "Cor between newt pheno and snake 0.976462242096197"
## [1] "pattern 1e-09_0.05_1e-10_0.158_0"
## [1] "Cor between average snake pheno and local cor 0.495645123018037"
## [1] "Cor between average newt pheno and local cor 0.683780730975857"
## [1] "Cor between average dif pheno and local cor 0.343077048380392"
## [1] "Cor between newt pheno and snake 0.863610636804903"
## [1] "pattern 1e-10_0.158_1e-09_0.05_0"
## [1] "Cor between average snake pheno and local cor -0.136365111793784"
## [1] "Cor between average newt pheno and local cor -0.00417878115632057"
## [1] "Cor between average dif pheno and local cor -0.403757876513809"
## [1] "Cor between newt pheno and snake 0.944867819529124"
## [1] "pattern 1e-10_0.5_1e-10_0.5_3"
## [1] "Cor between average snake pheno and local cor -0.18297822109376"
## [1] "Cor between average newt pheno and local cor -0.176414835554545"
## [1] "Cor between average dif pheno and local cor -0.106061513817461"
## [1] "Cor between newt pheno and snake 0.719531691596149"
## [1] "pattern 1e-09_0.158_1e-11_1.58_1"
## [1] "Cor between average snake pheno and local cor -0.159994859909367"
## [1] "Cor between average newt pheno and local cor -0.160427813829525"
## [1] "Cor between average dif pheno and local cor -0.00301056525820725"
## [1] "Cor between newt pheno and snake 0.821072242973623"
## [1] "pattern 1e-11_1.58_1e-08_0.05_1"
## [1] "Cor between average snake pheno and local cor 0.113070227081819"
## [1] "Cor between average newt pheno and local cor 0.156798285756727"
## [1] "Cor between average dif pheno and local cor -0.100138046198342"
## [1] "Cor between newt pheno and snake 0.847975459477899"
## [1] "pattern 1e-11_5_1e-09_0.5_0"
## [1] "Cor between average snake pheno and local cor -0.0155602209793962"
## [1] "Cor between average newt pheno and local cor 0.384223188661971"
## [1] "Cor between average dif pheno and local cor -0.164901108909709"
## [1] "Cor between newt pheno and snake -0.0409191734896981"
## [1] "pattern 1e-11_5_1e-08_0.158_2"
## [1] "Cor between average snake pheno and local cor 0.500230439052775"
## [1] "Cor between average newt pheno and local cor 0.36378156098346"
## [1] "Cor between average dif pheno and local cor 0.303737969267246"
## [1] "Cor between newt pheno and snake 0.634194505117167"
## [1] "pattern 1e-10_1.58_1e-09_0.5_3"
## [1] "Cor between average snake pheno and local cor 0.0539083740474133"
## [1] "Cor between average newt pheno and local cor 0.282863215735583"
## [1] "Cor between average dif pheno and local cor -0.202438537094216"
## [1] "Cor between newt pheno and snake 0.135596071275061"
## [1] 0.54326
## [1] 0.5474049
## [1] 0.5392002
## [1] 0.517508
## [1] 0.1655529
## [1] 0.2519408
## [1] 0.095895
## [1] 0.419727
## [1] 0.06514192
## [1] 0.3005453
## [1] 0.3867365
## [1] 0.4976603
## [1] 0.1991463
## [1] 0.08395684
## [1] 0.3101617
## [1] 0.5130782
## [1] -0.0288727
## [1] 0.1288655
## [1] -0.129634
## [1] 0.5223597
## [1] -0.0843053
## [1] -0.06692039
## [1] 0.1317734
## [1] 0.3610041
## [1] 0.2723125
## [1] 0.3114732
## [1] 0.3917085
## [1] 0.5393433
## [1] 0.2032271
## [1] -0.00932744
## [1] 0.1015125
## [1] 0.4523698